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1.  ZMTRODUCTZOII  AMD  SCOPE 


The  cumulative  probability  of  detection  of  Multi-Site  Daunage  fMSD) 
in  fuselage  lap-joints  of  aging  aircraft  was  assessed  in  a  previous 
study  [1].  It  showed  that  inspection  intervals  on  the  order  of  3000 
to  4000  flights  provide  cumulative  probabilities  of  detection  on 
the  order  of  0.98  (98  percent  detected)  or  better  if  inspections 
are  perfozmed  by  eddy  current.  Although  quite  a  few  assumptions 
were  necessary,  a  sensitivety  analysis  showed  the  effect  of  the 
assumptions  to  be  small. 

It  is  unlikely  for  MSD  to  occur  very  early  in  life,  because 
disbonding  of  the  adhesive  will  tedce  some  time.  Therefore,  and 
because  the  probability  of  detection  is  virtually  zero  in  the  early 
stages  of  MSD,  it  is  un-economical  to  adhere  to  the  prescribed 
inspection  interval  from  time  zero.  Hence  an  inspection  threshold 
was  instituted  by  stipulating  that  inspections  for  MSD  need  not 
start  until  after  the  accumulation  of  a  certain  number  of  flights. 
The  present  study  was  undertaken  to  provide  a  justification  for 
this  threshold. 

As  the  slotted  time  for  this  work  was  very  limited,  only  readily 
available  results  of  fatigue  tests  on  lap-joints  were  analyzed.  The 
test  data  were  used  primarily  to  obtain  the  scatter  (i.e. 
distribution  function)  for  the  start  of  cracking:  the  scatter 
determines  when  the  earliest  detectable  MSD  may  occur.  The  only 
true  service  experience  comes  from  the  Aloha  incident.  The  time  to 
failure  for  the  latter  case  was  used  as  the  anchor  point  for  the 
estimate  of  the  fatigue-life  cvirves  (S-N  curves)  for  service 
conditions.  As  this  approach  needed  only  a  few  reasonable 
assumptions,  it  is  considered  to  provide  a  realistic  assessment  of 
the  threshold. 

Another  procedure  to  determine  inspection  thresholds  is  based  upon 
crack  growth  analysis.  It  is  then  assumed  that  a  certain  crack  is 
already  present  in  the  new  structure,  and  the  time  it  tedces  for 
this  initial  crack  to  grow  to  a  "detectable”  size,  is  used  for  the 
inspection  threshold,  either  directly  or  with  a  safety  factor.  This 
procedure  has  found  some  acceptance  by  the  industry  and  by  FAA 
certification  offices.  The  initial  "crack"  is  not  a  real  crack  but 
rather  an  Equivalent  Initial  Flaw  (EIF) ,  where  the  qualifier 
"equivalent"  indicates  that  if  the  EIF  is  used  in  analysis,  it  will 
lead  to  acceptable  results.  Unfortunately,  the  results  depend 
strongly  upon  the  assumed  size  of  the  EIF,  while  the  commonly  used 
value  of  0.05  inch  cannot  very  well  be  justified.  Nevertheless, 
this  approach  was  used  in  the  present  work  as  well  to  obtain  a 
comparison  with  the  analysis  based  upon  fatigue  life.  The  results 
are  of  the  same  order  of  magnitude  as  those  based  upon  life 
analysis. 
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2.  KVhlLMMLM  UkTA 


2.1  DBSCSXPTZCai  OW  TOOK  BJkSB 

Readily  available  fatigue  test  data  for  joints  were  those  in 
References  2  through  8.  A  diligent  search  of  older  reports 
(1950-1965)  by  NACA  (NASA),  RAE,  NLR,  FFA,  ARL,  and  so  on,  is 
likely  to  provide  a  wealth  of  additional  data,  but  tiee 
constraints  did  not  permit  such  a  search  and  associated  data 
analysis.  Of  the  sources  readily  available  some  were  for  conditions 
and  configurations  not  immediately  relevant  to  the  problem  at  hand, 
but  the  l<u:ge  data  base  generated  by  Hartman  [  2  ]  is  very 
appropriate,  as  are  a  few  data  generated  recently  [6,7]. 

It  should  be  noted  that  in  all  of  the  following  the  stresses  are 
the  nominal  stresses  away  from  the  joint:  they  represent  the  hoop 
stress  in  a  fuselage.  This  simulated  hoop  stress  is  not  necessarily 
equal  to  that  due  to  presurization,  as  will  be  explained  later. 
Furthermore,  all  data  are  for  a  stress  ratio,  R,  of  0.05  to  0.1, 
while  fuselage  loading  is  essentially  at  R  «  0.  As  a  consequence, 
the  assumption  had  to  made  that  the  results  are  applicable  to  the 
slightly  different  fuselage  loading,  but  this  of  minor  importance 
because  the  effect  of  the  associated  difference  in  mean  stress  is 
small . 

Haurtman's  data  [2]  are  especially  useful,  because  they  were 
obtained  from  tests  on  adhesively  bonded  and  riveted  lap-joints  of 
a  configuration  (Figxire  1)  almost  identical  to  the  one  used  in  the 
fuselage  of  several  types  of  aircraft.  Hartman  performed  well  over 
400  tests.  The  lap-joints  specimens  were  bonded  with  a  cold-curing 
adhesive,  and  contained  3  rows  of  countersunk  fasteners.  Hartman 
investigated  the  effect  of  many  parameters,  the  most  important  of 
which  are:  (1)  Two  adhesives,  namely  AH106  (CIBA)  and  EC226  (3M) ; 
(2)  two  types  of  surface  treatment,  namely  chromic  acid  and 
sulphuric  acid  anodizing;  (3)  Three  testing  temperatures,  namely 
-55C  (stratosphere) ,  20C  and  50C  (aircraft  taking  of  after  exposure 
to  summer  sun  at  airport);  (4)  in  one  series  of  tests  the  adhesive 
was  intentionally  maltreated  by  exposure  to  100%  RH  at70C  for  4 
weeks;  (5)  three  loading  frequencies,  n<‘>mely  6,  240  and  3000  cpm; 
(6)  several  different  stress  levels. 

Part  of  the  specimens  failed  at  the  adhesive  fillet  or  elsewhere 
away  from  the  fastener  holes.  These  data,  regardless  of  test 
parameter,  are  shown  in  Figure  2.  As  MSD  occurs  at  the  fastener 
holes,  these  data  were  excluded  from  the  following  analysis.  Well 
over  200  of  Hartman's  specimens  failed  at  the  fastener  holes.  They 
are  shown  in  Figure  3,  again  regardless  of  test  parameters.  Hartman 
did  a  limited  number  of  tests  on  specimens  that  were  riveted  only 
(as  opposed  to  riveted  and  bonded) .  The  data  of  these  tests  are 
shown  in  Figure  4.  Also  shown  in  Figure  4  are  the  data  points 
obtained  at  low  frequency  loading  (which  is  the  most  relevant  for 
longitudinal  fuselage  joints)  for  riveted  plus  bonded  joints  which 
failed  at  the  holes. 
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The  data  obtained  by  Mayville  [6,7]  deserve  attention.  Instead  of 
the  basic  lap-joints  as  used  by  Hartman,  Mayville  employed  the 
specimen  shown  in  Figure  5.  Two  short  stiffeners  were  attached  to 
the  edges  of  the  specimens  to  simulate  the  crack  arrester  straps 
(tear  straps)  found  in  some  types  of  aircraft  experiencing  MSD.  A 
finite  element  analysis  [9]  had  shown  the  stresses  at  the 
lap-slplce  to  be  higher  midway  between  the  straps.  Strain  gage 
measurements  showed  that  the  stress  distribution  in  the  test  panels 
with  the  simulated  tear  straps  was  nearly  identical  to  the  one 
obtained  from  the  finite  element  analysis.  This  is  shown 
convincingly  in  Figure  5. 

The  test  data  from  these  panels,  although  limited,  may  be  more 
representative  of  the  situation  in  some  aircraft,  but  it  should  be 
noted  that  the  joints  were  riveted  only  (not  riveted  and  bonded) , 
and  that  they  were  obtained  at  high  frequency  loading.  A  total  of 
10  data  points  was  available;  in  €  of  the  specimens  the  fasteners 
had  driven  heads  of  nominally  0.24  inch  diameter,  while  in  the 
others  this  dimension  was  0.25  inch.  The  data  points  for  all  are 
shown  in  Figure  6.  Those  for  the  small  driven  heads  fall  below 
Hartman's  scatter  band;  the  others  fall  inside  the  scatterband. 
This  indicates  that  the  non-uniform  stress  distribution  in  the 
specimens  (Figure  5)  has  but  little  effect.  Figxire  6  shows  average 
lives  for  various  conditions  investigated  by  Hartman.  In  addition 
a  line  is  shown  labeled  "Grover"  [8].  Although  Grover  [8]  does  not 
identify  his  sources,  nor  the  details  of  the  specimen  used,  at 
least  the  curve  falls  within  the  scatterband  of  Hartman's  data  in 
the  regime  of  stresses  relevant  to  the  problem  of  longitudinal 
fuselage  lap-joints. 


2.2  SCATTER  AMD  DISTRIBUTION  FUNCTION 

As  is  normal  in  the  case  of  fatigue  (especially  at  lower  stresses) , 
the  scatter  in  Figures  3  and  4  is  large.  Admittedly,  this  scatter 
is  somewhat  exaggerated  because  all  data  for  all  conditions  are 
plotted  in  the  same  graph;  the  justification  for  this  will  follow. 
Despite  the  large  number  of  tests,  the  number  of  data  points  per 
stress  level  and  per  case  is  still  very  small  and  would  yield 
distribution  functions  containing  much  uncertainty.  This  problem 
was  solved  by  the  use  of  the  fatigue  damage  at  failure  according  to 
Miner's  rule,  n/N,  where  n  is  the  actual  number  of  cycles  to 
failure,  and  N  the  average  life  at  the  same  stress  level.  The 
advantage  of  this  is  that  data  for  all  stress  levels  are 
consolidated.  To  obtain  the  distribution  function,  only  those  data 
were  considered  that  pertain  to  failure  at  the  rivet  holes  (Figures 
3,4,6),  as  these  are  the  most  relevant  to  MSD. 

This  leaves  about  200  data  points  for  which  the  distributions  are 
shown  in  Figure  7.  Scatter  for  the  individual  cases  is  very  similar 
to  the  total  scatter  in  Figure  8  indicated  by  the  fact  that  they 
all  go  from  n/N  -  0.2  to  2  (plus).  Besides,  the  distribution 
functions  are  the  same,  as  all  data  follow  the  same  curve 
reasonably  well.  In  the  case  of  actual  aircraft  almost  all 
conditions  investigated  by  Hartman  will  appear  part  of  the  time. 
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Therefore  all  data  were  considered  to  belong  to  one  population 
covering  all  service  conditions  for  aircraft.  This  leads  to  the 
distribution  function  sho%m  in  Figure  8.  Not  only  is  this  general 
distribution  fvmction  the  same  as  the  one  for  the  individual  cases 
in  Figure  1 ,  it  turns  out  that  the  data  fit  the  Weibul  distribution 
very  well. 

Mi''er's  rule  predicts  crack  inititation  if  Sum  n/N  »  1.  The  data 
show  that  due  the  scatter,  failure  occurs  at  values  between  0.15 
and  2.  The  Weibul  parameters  in  Figure  8  were  more  or  less 
confirmed  by  other  data  [3,4,5].  The  distribution  function  is: 


Pf  «  1  -  exp(-  (n/N  -  Do) /(A  -  Do) x) 

where  Pf  is  the  fraction  (percentage)  failed,  or  probability  of 
failure,  while  Do,  A,  and  X  are  the  shape  parameters  the  values  of 
which  are  shown  in  Figures  7  and  8.  The  average  (Pf  »  0.5)  is 
indeed  at  n/N  »  1. 
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3.  DZSTRZBUTiVOM  07  ZVZTZATZOII  LZTB  ZV  PRKCTZCB 


3.1  BBTZIOkTB  07  BPPLZCABLB  8-lf  CURVE 

Although  the  above  provides  the  distribution  of  n/N,  the  actual 
average  life,  N,  under  aircraft  service  conditions  is  not  known. 
However,  one  important  actual  service  data-point  is  available  from 
the  Aloha  incident:  the  failure  occured  at  89600  flights. 
Accounting  for  the  fact  that  the  average  time  for  crack  growth  to 
failure  is  25,000  flights  the  initiation  time  is  estimated  at 
64,000  flights. 

The  average  of  25,000  flight  cycles  for  crack  propagation  follows 
from  previous  work  [10]  in  which  crack  growth  was  calculated  for  a 
variety  of  statistical  parameters  (assigned  by  a  Monto-Carlo 
technique)  and  for  a  variety  of  circumstances.  An  example  of  these 
computations  is  shown  in  Figure  9.  The  number  is  confirmed  by  data 
obtained  by  Mayville  [6,7]  shown  in  Figure  10.  Although  the  crack 
growth  curves  in  Figure  10  start  at  different  cycle  nximbers,  they 
are  well-nigh  parallel.  The  total  growth  from  initiation  to  failure 
covers  about  20,000  to  25,000  cycles,  if  initiation  is  considered 
as  the  appearance  of  a  crack  of  about  0.02  inch  size.  From  a 
practical  point  of  view  this  is  a  good  definition,  because  at  0.02 
inch  the  crack  will  just  emerge  from  under  the  fastener  head  and 
become  inspectable. 

The  remaining  question  is  then  where  the  Aloha  case  falls  in  the 
scatter  band,  i.e.  whether  the  case  is  an  extreme.  One  might  argue 
that  it  is,  because  Aloha  operates  in  adverse  conditions  (low 
altitude  flights  over  sea) .  But  USD  was  detected  in  many  other 
aircraft  as  well,  indicating  that  Aloha  was  not  an  extreme.  Three 
cases  were  considered,  namely: 

1.  Aloha  is  average  (or  50  percentile;  conservative). 

2  Aloha  belongs  to  the  lower  25  percentile  . 

3.  Aloha  belongs  to  lower  15  percentile  (optimistic;  it  is 
an  extreme) . 

For  each  of  these  percentiles  the  value  of  n/N  can  be  read  from 
Figure  8,  from  which  the  average  life  N  can  be  found.  Making  the 
reasonable  assumption  that  all  S-N  curves  have  the  same  slope,  the 
life-to- initiation  curves  for  the  above  3  conditions,  and  for  the 
stress  levels  encountered  in  service  can  be  determined  to  be  as 
shown  in  Figure  11,  and  as  explained  below. 

For  the  relevant  stress  ranges  in  the  regime  of  10  to  20  ksi  the  5- 
N  curve  can  be  represented  by  a  straight  line  on  semi-logarithmic 
scales,  and  hence  by  the  equation: 


S  =  A  -  q  In  N 
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which  can  be  inverted  to: 

N  -  exp  ((A  -  S)/q) 


where  N  is  the  life  to  crack  initiation,  S  the  stress  range,  and  A 
and  g  are  parameters.  That  all  lines  should  be  parallel  (same 
slope)  in  the  region  of  interest  can  be  ascertained  from  Figures  3, 
4,  and  6;  this  results  in  a  fixed  value  for  g  of  2.74.  (Although 
the  data  in  the  previous  figures  were  plotted,  in  accordance  with 
engineering  practice,  on  10>log  scale,  the  natural  logarithm  was 
used  in  the  equations.)  Hence,  the  difference  in  the  lines  shotm  in 
Figure  11  is  characterized  by  different  values  of  A. 

As  flights  are  of  different  length  the  cruising  altitude  (euid  hence 
the  differential  pressure)  varies  from  flight  to  flight.  This  means 
that  1-he  loading  is  not  of  constant  amplitude.  Since  the  5-N  curves 
were  to  be  "anchored"  on  the  Aloha  case,  the  Aloha  experience 
service  experience  is  detailed  below: 


Altitude  Delta-p 
(ft)  (psi) 


%  of  flights  Pressure  relative  to  highest 


>18000 

7.5 

23 

16000 

6.7 

5 

<14000 

6.1 

72 

1.00 

0.89 

0.81 


According  to  the  finite  element  analysis  [9]  the  membrane  stress  at 
mid-bay  is  14.9  ksi,  but  this  number  is  for  a  skin  thickness  of 
0.04  inch  and  a  differential  pressure  of  8.5  psi.  In  the  actual 
aircraft  the  maximum  pressure  was  7.5  psi  (above  table),  bringing 
the  maximim  operating  stress  to  14.9  *  7. 5/ 8. 5  «  13.15  ksi.  Due  to 
the  lower  skin  thickness  another  10  percent  must  be  added,  which 
provides  a  maximum  stress  of  14.5  ksi. 


If  the  Aloha  case  was  the  average  (50  %)  crack  inititiation,  as 
defined  above,  occurred  at  n/H  »  1,  at  64,000  cycles  as  explained. 
The  table  shows  that  the  highest  stress  of  14.5  ksi  occurred  in  23 
percent  of  the  flights,  l.e.  in  0.23  *  64,000  «  14,720  flights. 

The  membrane  stress  in  the  other  flights  at  differential  pressures 
of  6.7  and  6.1  psi  are  obtained  as  14.5*6.7/7.5  «  13  ksi,  and  14.5 
*  6. 1/7.5  s  11.8  ksi,  respectively.  The  associated  cycle  numbers 
are  0.05  *  64,000  3,200,  and  0.72  *  64,000  «  46,000, 
respectively . 

With  this  information  a  system  of  4  equations  with  4  unknowns  can 
be  solved  to  obtain  A.  The  result  is  A  43,  so  that  the  equation 
for  the  life  becomes: 


N  *  exp((43  -  S)/2.74) 
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The  neabrane  stress  in  the  other  flights  at  differential  pressures 
of  6.7  and  6.1  psi  are  obtained  as  14.5*6.7/7.5  «  13  ksi,  and  14.5 
*  6. 1/7. 5  ■>  11.8  ksi,  respectively.  The  associated  cycle  numbers 
are  0.05  *  64,000  *  3,200,  and  0.72  *  64,000  »  46,000, 
respectively. 

With  this  information  a  system  of  4  equations  with  4  unknowns  can 
be  solved  to  obtain  A.  The  result  is  A  «  43,  so  that  the  equation 
for  the  life  becomes: 

N  =  exp((43  -  S)/2.74) 

Inverting  this  procedure  permits  calculation  of  sum  n/N  for  this 
case  as  follows: 


(ksi) 

n 

N  (for  A  ®  43) 

n/N 

14.5 

14,720 

32,900 

0.447 

13.0 

3,200 

56,900 

0.056 

11.8 

46,000 

88,200 

0.522 

Sum  n/M 

1.025 

;  Aloha 

is  an  average 

case,  the  S-N  curve 

with  A 

shown  in  Figure  11,  is  the  applicable  curve  for  the  average  life  to 
initiation.  For  the  other  cases  defined  above  a  similar  procedure 
leads  to  the  average  S-N  curve.  However,  it  can  be  seen  immediately 
what  the  results  will  be,  because  all  calculations  are  based  upon 
proportionality.  Figure  6  shows  that  for  a  probability  of 
initiation  of  25  percent  the  value  of  n/N  =  0.70,  while  for  a 
probability  of  15  percent  it  is  0.52.  Thus  it  can  be  deduced 
immediately  that  the  associated  lives  for  all  stresses  are  higher 
by  a  factor  of  1/0.70  =  1.43,  and  1/0.52  =  1.92  respectively, 
leading  to  values  for  A  of  44  and  44.8  respectively;  these  are 
represented  by  the  other  lines  in  Figure  11. 


Although  superfluous,  it  should  be  pointed  out  that  for  the 
following  computations  the  lines  in  Figure  11  should  be 
interpreted  as  averages.  For  example,  if  the  line  for  A  »  43  is  the 
average,  the  Aloha  case  will  be  average  (as  signified  by  the  Aloha 
case  falling  on  this  line) ;  but  if  the  line  for  A  =  44  is  the 
average,  the  Aloha  case  falls  left  of  the  line,  as  shown  in  Figure 
11.  The  curves  in  Figiire  11  pertain  to  the  average  life  N.  To 
this  the  scatter  (distribution  function)  of  Figure  8  was  applied. 


3.2.  CALCULATION  OF  INITIATION  LIFE  AND  SCATTER 

The  S-N  curve (s)  and  distribution  function  now  being  known,  the 
life  to  initiation  can  be  calculated.  To  this  end  a  small  computer 
program  was  developed.  The  program  employs  the  membrane  stresses 


across  a  bay  as  calculated  by  the  finite  elenent  analysis  [9]  for 
a  pressure  differential  of  8.5  psl  (Figure  5),  as  shown  below: 


11.1  ksl  (at  frame  and  straps)  at  10%  of  holes 

12.6  ksl  at  20%  of  holes 

13.8  ksl  at  20%  of  holes 

14.5  ksl  at  20%  of  holes 

14.8  ksl  at  20%  of  holes 

14.9  ksl  (midbay)  at  10%  of  holes 

The  nominal  membrane  stress  was  used  in  the  analysis,  because  all 
data  previously  discussed  are  In  terms  of  this  nominal  stress.  The 
nominal  stress  for  any  pressure  differential  the  can  be  obtained 
by  multiplying  the  above  stresses  by  p/8.5. 

Many  airlines  service  longer  routes  than  Aloha,  so  that  their 
maximum  pressxire  (altitude)  will  be  higher.  It  seems  reasonable 
however,  lacking  data  from  other  airlines,  to  assume  that  the  last 
two  columns  in  the  Aloha  usage  table  provide  an  appropriate 
estimate  of  "general”  usage  relative  to  maximum  pressure 
differential.  Therefore,  the  calculations  were  based  upon  a  usage 
spectrum  in  accordance  with  the  one  shown  for  Aloha,  relative  to  a 
maximum  operating  pressure  of  1,  and  Miner's  rule  for  the 
accumulation  of  damage  at  the  three  stress  ranges  (differential 
pressures) .  The  maximum  operating  pressure  differential  is  then 
the  only  variable,  because  the  other  pressures  are  given  relative 
to  the  maximum.  It  permits  computation  of  the  membrane  stresses  at 
the  holes  according  to  the  conversion  rule  discussed,  and 
subsequently,  assessment  of  the  damage,  Siun  n/M,  according  to  the 
generalized  usage  spectrum. 

Apart  from  being  directly  dependent  upon  fuselage  pressurization, 
the  membrane  stress  depends  upon  aerodynamic  pressure.  The 
aerodynamic  pressure  varies  from  point-to-point  due  to  the  fuselage 
shape,  but  should  be  assessed  at  about  0.5  psl  on  average  [12]  . 
This  will  increase  the  pressure  differential  locally,  and  hence, 
the  membrane  stresses  will  be  higher  than  those  following  from 
fuselage  pressurization.  This  may  explain  why  MSD  appears  to  be 
more  prevalent  at  joints  in  the  fuselage  crown  just  behind  the 
flight  deck.  Another  area  where  this  effect  is  significant  is  the 
vicinity  of  the  wing-fuselage  connection. 

If  the  pressure  differential  by  fuselage  pressurization  is  for 
example  8  psi,  the  actual  local  pressure  differential  would  be  on 
the  order  of  8.5  psi,  and  possibly  higher  close  to  the  wing- 
fuselage  connection.  The  above  was  invoked  in  the  small  computer 
program  already  mentioned. 

Accounting  for  all  the  effects  discussed  above  calculations  were 
made  of  the  per  cent  failed  (holes  cracked) .  The  computations  were 
cintinued  until  10  percent  of  the  holes  were  cracked,  for  the 
cases  discussed  and  for  different  pressures.  The  results  of  the 
computations  are  shown  in  Figures  12  through  14  for  various 
pressure  differentials  (as  adjusted  for  usage  spectrum  and  local 
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aerodynamic  pressure) ,  and  for  cases  where  the  Aloha  case  is  the 
lower  15,  25  and  50  percentile.  For  a  particular  maximum  pressure 
differential  of  8  psi  (accounting  for  an  aerodynamic  differential 
of  0.5  psi),  the  results  are  shown  in  Figure  15.  Similar  cross 
plots  can  be  made  for  other  pressures  by  using  the  data  in  Figures 
12  through  14.  The  results  are  trivial  in  a  way,  because  they 
merely  restate  the  lower  end  of  the  distribution  function  for  a 
particular  set  of  circumstances. 
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4.  COMSZDERlkTZOMS  FOR  THRBSBOU) 


4.1.  DBFZHZTZON  OF  THRESHOLD 

A  more  specific  definition  of  threshold  is  now  needed.  In  essence 
the  threshold  is  the  tine  at  which  inspections  must  be  started. 
There  is  no  use  for  inspections  as  long  as  the  cracks  cannot  be 
found,  i.e.  if  they  are  so  small  that  the  probability  of  detection 
is  virtually  zero.  From  execution  of  the  computer  program  for  crack 
detection  [10],  it  is  known  that  the  largest  crack  is  on  the  orderr 
of  0.06  inch  when  2  percent  of  the  holes  is  cracked  and  0.12  inch 
when  5  percent  is  cracked;  all  other  cracks  are  smaller.  In  view  of 
that  the  probability  of  detection  is  virtually  zero.  Although  the 
results  can  vary  considerably  in  different  computer  runs  (the  code 
simulates  statistical  variability  by  means  of  a  Monte-Carlo 
technique  [10]),  the  above  numbers  are  considered  reasonably 
conservative . 

Since  the  following  argiiments  will  be  based  upon  the  results  of  a 
computer  program  developed  previously  [10],  it  seems  appropriate  to 
demonstrate  that  the  computations  provide  results  in  concert  with 
actual  observations  of  MSD.  Figure  16  shows  an  exsuaple  of  computed 
crack  sizes  at  100  holes  at  the  time  of  failure,  while  Figxire  17 
shows  the  computed  cumulative  probability  of  detection  of  these 
cracks.  Figure  16  should  be  compared  with  Figures  18  and  19.  The 
latter  two  figures  show  the  MSD  crack  sizes  as  observed  by  Mayville 
[6,7]  in  test  panels  (Figxire  18),  and  as  detected  in  an  aircraft 
(Fig\ire  19) .  Apparently,  the  computations  produce  a  "true-to-life” 
picture. 

The  computer  program  was  therefore  used  to  produce  the  MSD  crack 
sizes  for  the  situation  where  the  largest  crack  is  0.06  inch,  and 
0.1  inch  respectively.  Only  two  examples  will  be  provided.  Figure 
20  shows  a  case  where  the  fraction  of  holes  cracked  is  0.09  (9 
percent  cracked) ,  while  Figure  21  shows  the  growth  curve  of  the 
largest  crack  as  well  as  the  number  of  cracks  as  a  function  of  the 
number  of  (flight)  cycles.  As  already  mentioned,  the  results  vary 
considerably  from  run  to  run.  A  rather  extreme  case  in  which 
already  39  percent  of  the  holes  are  cracked  -  the  largest  crack 
still  being  0.1  inch  -  is  shown  in  Figure  22. 

More  important  than  Figures  20  and  22  would  be  the  figures  showing 
the  cumulative  probability  of  detection  (in  the  manner  shown  in 
Figure  17) .  However,  these  figures  would  exhibit  the  "scales"  only, 
because  the  cumulative  probability  of  detection  was  virtually  zero 
in  all  cases  and,  therefore,  would  not  show  in  the  graphs.  These 
calculations  were  done  for  an  inspection  interval  of  3000  flights. 
Had  the  interval  been  selected  larger  than  this,  the  probability  of 
detection  would  have  been  less  (if  less  than  zero  were  possible) . 
Naturally,  shorter  intervals  would  show  higher  cumulative 
probability  of  detection,  but  that  is  of  academic  interest  only, 
because  the  interval  is  not  less  than  3000  flights. 
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4.2  RB8ULTIM6  FLIGHT  NUMBERS  TILL  THRESHOLD 


The  results  indicate  that  the  definition  of  threshold  nay  well  be 
the  tine  at  which  5  percent  of  the  holes  are  cracked,  because  the 
cunulative  probabilty  of  detection  would  still  be  virtually  zero. 
However,  it  is  not  the  charter  of  the  author  of  this  report  to 
suggest  what  the  threshold  should  be.  This  report  nerely  provides 
infomation  upon  which  the  authorities  can  base  a  decision. 

For  this  reason  thresholds  defined  by  1,  2  and  5  percent  cracked 
were  considered.  The  nunber  of  flights  for  reaching  these  can  be 
obtained  fron  the  basic  results  provided  in  Figures  12  through  14. 
To  facilitate  interpretation  additional  cross  plots  were  nade  for 
the  percentages  nentioned  above.  These  are  shown  in  Figures  23 
through  26. 

The  following  nay  serve  as  an  exanple  of  how  these  figures  should 
be  interpreted.  If  one  is  willing  to  assune  that  the  Aloha  case 
belongs  to  the  lower  25  percentile,  that  the  naxinun  operating 
pressure  is  8  psi,  and  that  5  percent  cracked  is  a  good  definition 
of  threshold,  the  resulting  threshold  would  be  34,000  flights 
(Figures  25  and  26) .  With  very  conservative  assunptions  (Aloha  is 
average,  2  percent  cracked,  naxinun  operating  pressure  of  9  psi), 
the  threshold  would  be  10,000  flights  (Figure  24). 


4.3  THE  CRACK  GROWTH  APPROACH 

Inspection  thresholds  are  sonetines  determined  by  neans  of  crack 
growth  analysis.  In  that  case  a  certain  initial  crack,  denoted  as 
the  Equivalent  Initial  Flaw  (EIF) ,  is  assuned  present  in  the  new 
structure.  The  tine  (nunber  of  cycles)  it  takes  for  this  EIF  to 
grow  to  a  detectable  size  is  used  as  the  basis  for  the  inspection 
threshold.  The  EIF  is  often  taken  as  0.05  inch;  this  size  is  based 
upon  a  sonewhat  arbitrary  EIF  detemined  by  the  USAF,  as  explained 
below. 

The  USAF  [11,12]  exanined  the  cracks  that  developed  during  a  full- 
scale  fatigue  test  on  an  F-4  wing.  Of  the  2000  holes  present  a 
total  of  119  had  developed  cracks.  As  the  load  history  and  stresses 
in  the  test  were  known,  calculations  could  be  nade  of  the  growth  of 
those  119  cracks.  The  calculations  were  adjusted  to  natch  the 
final  crack  sizes  observed  at  the  end  of  the  test.  It  turned  out 
that,  in  order  for  these  cracks  to  have  developed  to  the  size  at 
the  end  of  the  test,  the  calculations  would  have  to  assvine  that  a 
certain  initial  was  already  present  in  the  new  structure.  This 
initial  flaw  was  clearly  an  equivalent  flaw,  which  would  nake  the 
conputations  conpatible  with  the  final  crack  observed.  The  EIF 
derived  fron  these  conputations  was  on  the  order  of  a  few  nils. 
Taking  the  distribution  of  the  calculated  sizes  of  the  EIF  for  the 
119  holes  cracked  (the  1881  uncracked  holes  which  would  have 
provided  an  EIF  of  zero  were  ignored) ,  an  estinate  was  nade  of  the 
EIF  needed  for  extrene  probability  of  occurrence.  While  a  nomal 
distribution  and  a  Weibul  distribution  would  have  yielded  a  nuch 
snaller  EIF,  a  Johnson  distribution  was  assuned,  which  exaggerates 


extreme  values.  On  top  of  that  the  1881  that  did  not  develop  cracks 
(EIF  »  0)  were  ignored.  These  assumptions  led  to  an  EIF  of  0.05 
inch,  which  -in  view  of  the  above-  is  a  rather  arbitrary  size. 

Subsequently,  many  fatigue  tests  were  performed  [12,13]  on 
specimens  with  holes  to  substantiate  the  0.05  inch  EIF. 
Invariably,  these  tests  showed  that  the  EIF  is  on  the  order  of 
0.001  to  0.002  inch.  Be  that  as  it  may,  if  the  USAF  achieves  safe 
aircraft  by  assuming  an  EIF  of  0.05  inch,  the  assumption  cannot  be 
argued  with.  Unfortunately  however,  the  0.05  inch  EIF  is  often 
considered  the  final  answer,  and  used  in  non-military  damage 
tolerance  analysis  as  if  it  has  a  sound  basis. 

Despite  the  above,  a  crack  growth  analysis  based  upon  an  EIF  was 
used  in  the  present  work  to  obtain  an  estimate  of  the  threshold. 
There  is  one  obvious  problem  however.  As  shown  in  the  previous 
section,  the  crack  size  is  already  0.05  inch  when  2  percent  of  the 
holes  are  cracked.  Hence,  a  crack  growth  analysis  starting  with  an 
EIF  of  0.05  would  yield  no  life,  and  would  lead  to  an  inspection 
threshold  of  zero,  if  the  definition  of  threshold  were  2  percent 
cracked.  This  can  be  easily  ascertained  from  previous  figiires: 
growth  from  a  crack  size  of  0.05  inch  to  the  threshold  crack  size 
of  0.1  inch  covers  only  a  few  thousand  cycles.  Therefore  the 
analysis  was  based  upon  a  smaller  EIF  (Figure  27),  and  the  results 
were  inverted  to  show  the  life  for  other  values  of  the  EIF  (Figure 
28) . 

A  comparison  of  Figures  27  and  28  with  Figures  9  and  10  shows  that 
the  computed  cycle  numbers  are  realistic  and  believable  (note  that 
the  curves  in  Figure  10  are  for  a  pressure  differential  of  8.5  psi, 
and  should  be  compared  with  those  for  a  pressure  differential  of  8 
psi  in  Figures  27  and  28,  because  the  addition  of  0.5  psi 
aerodynamic  pressure  brings  the  differential  pressure  at  8.5  psi). 

It  is  obvious  from  Figure  28  that  the  assumption  regarding  the  size 
of  the  EIF  is  crucial  for  the  result.  It  is  not  the  charter  of  the 
author  of  this  report  to  make  recommendations  or  suggestions. 
Therefore,  the  results  are  presented  "as  is”,  so  that  someone  in 
authority  can  draw  conclusions  after  having  decided  upon  a  suitable 
EIF  size. 
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5.  CQHCLU8ZOM 


The  threshold  for  the  start  of  inspections  for  MSD  in  longtudinal 
fuselage  lap  joints  will  be  on  the  order  of  15,000  -  30  000 
flights,  depending  upon  the  definiton  of  threshold.  It  was  shown 
that,  even  if  5  percent  of  the  holes  are  cracked,  the  probability 
of  detection  of  the  MSD  is  virtually  zero.  Therefore,  the 
definition  of  threshold  may  well  be  the  number  of  flights  at  which 
5  percent  of  the  holes  is  cracked,  in  idiich  case  the  larger  of  the 
numbers  quoted  applies.  Provided  one  makes  the  appropriate 
assumption  for  the  size  of  the  equivalent  initial  flaw  present  in 
the  new  structiire,  a  crack  growth  analysis  leads  to  approximately 
the  same  conclusion. 
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FIGURE  1.  TESTS  BY  HARTHAN 


FIGURE  2.  ALL  DATA  FOR  FAILURES  IN  ADHESIVE;  NOT  AT  HOLES  -  ALL  CONDITIONS 
(HARTMAN  DATA  FOR  RIVETED  PLUS  BONDED  LAP-JOINTS) 


(IS>I)  39NVH  SS3«iS 
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FIGURE  3.  ALL  DATA  FOR  FAILURES  AT  HOLES  (HARTMAN  DATA  FOR  RIVETED  AND 

BONDED  LAP-JOINTS;  HARTMAN  DATA) 


LOW  FREQUENCY  AND  HOLE-FAILURES 
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FIGURE  4.  HARTMAN  DATA  FOR  LAP-JOINTS  RIVETED  PLUS  BONDED  AND  RIVETED  ONLY 

(CONSOLIDATION  OP  ALL  HARTMAN  DATA) 
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FIGURE  5.  TESTS  ON  SPECIMEN  WITH  SIMULATED  TEAR  STRAPS  t6»7] 
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FIGURE  6.  COMPARISON  OF  MEAN  S-N  CURVES 


CURVE;  Do  -  0.15:  LAMBDA  »  1.25;  ALPHA  -  1.73 

>  -55C:  NOT  AGED:  HF 
0  -55C:  AGED  :  HF 
X  20C:  NOT  AGED:  HF 
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FIGURE  7.  HOLE-FAILURES  AT  HIGH  AND  LOW  FREQUENCY  -  ADHESIVE  AND  NON-AGED 
(HARTMAN  DATA  FOR  RIVETED  PLUS  BONDED  LAP-JOINTS) 


FIGURE  9.  GROWTH  OP  LARGEST  CRACK  IN  10  RUNS 
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FIGURE  10.  MEASURED  CRACK  GROffTH  CURVES  [6,7] 
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FIGURE  11.  S-N  CURVES  USED  IN  ANALYSIS 


MAXIMUM  DIFFERENTIAL  PRESSURE 
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.SE  WHERE  ALOHA  BELONGS  TO  LOWER  15  PERCENTILE 


MAXIMUM  DIFFERENTIAL  PRESSURE  IN  psi 
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FIGURE  13.  CASE  IN  WHICH  ALOHA  BELONGS  TO  LOWER  25  PERCENTILE 


MAXIMUM  DIFFERENTIAL  PRESSURE 
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FIGURE  14.  CONSERVATIVE  CASE  WHERE  ALOHA  IS  AVERAGE 
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FIGURE  15.  THREE  CASES  FOR  DIFFERENTIAL  PRESSURE  OF  8  psl 
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FIGURE  16.  FINAL  CRACK  SIZES  AT  ALL  LOCATIONS  -  SAME  CASE  AS  IN  FOLLOWING 
FIGURE  -  SOLID  LINE:  LEFT  CRACK;  DASH-DOT:  RIGHT  CRACK 
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FIGURE  17.  CUMULATIVE  PROBABILITY  OF  DETECTION  PER  LOCATION  -  SOLID  LINE:  EDI 
DASH-DOT  LINE:  VISUAL  INSPECTION  INTERVAL  3000  CYCLES;  TYPICAL  CASE 
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FIGURE  18.  MEASURED  CRACK  SIZES  AT  FAILURE  [6,7] 


FIGURE  19.  SIZES  OF  DETECTED  NSD  IN  AIRCRAFT  [6,7] 


FIGURE  20.  SIZE  DISTRIBUTION  WHEN  LARGEST  CRACK  IS  0.1  INCH 


FIGURE  21. NUMBER  OF  CRACKS  AND  SIZE  OF  LARGEST  CRACK  -  EXAMPLE  OF  STATISTICALLY 
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FIGURE  22.  SIZE  DISTRIBUTION  WHEN  LARGEST  CRACK  IS  0.1  INCH 
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FIGURE  23.  SUMMARY  OF  RESULTS  -  LIFE  UNTIL  1  PERCENT  CRACKED 
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MAX.  PBESSUHE  DIFF.  (psi) 
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FIGURE  25.  SUMMARY  OF  RESULTS  -  LIFE  UNTIL  5  PERCENT  CRACKED 


FIGURE  26.  SUMMARY  OF  RESULTS  -  CASE  WHERE  AIX)HA  BEIX)NGS  TO  IX)WER  25  PERCEHTILB 
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STARTING  AT  0.005  INCH 
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FIGURE  28.  LIFE  UNTIL  0.12  INCH  CRACK  AS  A  FUNCTION  OF  ASSUMED  EQUIVi 

INITIAL  FLAW  (EIF)  SIZE 


